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Abstract; We analyze an ensemble in which energy (£), temperature (T) and multiplicity (A'^) can all fluctuate and with the 
help of nonextensive statistics we propose a relation connecting all fluctuating variables. It generalizes Lindhard's thermodynamic 
uncertainty relations known in literature. 
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1 Introduction 

A long time ago, it was suggested that the thermodynamical quantities, temper- 
ature T and energy could be regarded as being complementary in the same 
way as are position and momentum in quantum mechanics The reason is 
that the only way to attribute a definite temperature to a physical system is 
by bringing it into thermal contact and equilibrium with another very large 
system acting as a heat bath. In this case, however, the system will freely ex- 
change energy with the heat bath, and one loses the possibility of controlling 
its energy. On the other hand, in order to make sure that the system has a 
definite energy, one should isolate it from its environment. But then there is 
no way to determine its temperature. Dimensional analysis already leads to the 
conjecture that this relation would take the form AU A(3 > k, where f5 = 1/T 
and k is Boltzmann's constant. The isolation [U definite) and contact with a 
heat bath (T definite) then represent two extremal cases of such complementar- 
ity. This idea has so far not received much recognition in the literature because 
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its validity and foundations are under discussion (see [2] for review and [3] for 
comments). The main point is the exact meaning of the increments A. Indeed, 
all versions of this uncertainty relation proposed so far employ different theo- 
retical frameworks and give different interpretations of the uncertainty A/3, in 
most cases by using concepts from theories of statistical inference. 

In this paper we shall treat these increments as a measure of fluctuation of the 
corresponding physical quantities and analyze an ensemble in which energy (t/), 
temperature (T) and multiplicity {N), can all fluctuate. In this way we gener- 
alize the relation between fluctuations of U and T derived in thermodynamics 
[3] expressed by their relative variances Var{x) / {x^ = as: 



Eq. ([T]) is an attempt at describing a small system remaining in thermal contact 
with a heat bath of varying size. It represents the kind of uncertainty relation 
mentioned before, namely that the standard deviation of one variable can be 
made small only at the expense of increasing the corresponding standard devia- 
tion of the conjugate variable [2113] • This relation is supposed to be valid all the 
way from the canonical ensemble, for which Var{T) = and Var{U) = 1/ (N), 
up to the microcanonical ensemble for which Var{T) = 1/{N) and Var{U) = 0. 
It expresses the complementarity between the temperature and energy and the 
canonical and microcanonical description of the system. The generalization pro- 
posed here extends relation (jlD, also including in it the possible fluctuations of 
multiplicity and its connections with fluctuations of temperature; this is done 
using nonextensive Tsallis statistics [3] which incorporates such fluctuations in 
a natural way PI^. We shall demonstrate that such an approach results in a 
characteristic relation connecting all fluctuating variables. 

2 Fluctuations in statistical mechanics 

2.1 Single-variable fluctuations 

Suppose that in a process one has independently produced secondaries with 
energies {£'i=i,...,iv}, each distributed according to Boltzmann distribution char- 
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acterized by a temperature T 0, 

5, (E,) = ^ • exp • E,) where p-^ = T. (2) 
The corresponding joint probabihty distribution in this case is given by 

( N \ N 

9 {{E^^i....,n]) = /^^ exp -/^ E ^0 = ^^P(-^ U=Y.E,. (3) 

V i^l J i=l 

When two of three variables {U,T,N) are fixed, one has three possible situa- 
tions: 

• T = 1/13 and are fixed and the energy U (because {£'i=i,...,iv} are indepen- 
dent) can fluctuate. Using characteristic functions or sequentially performing 
integration of the joint distribution (I3j) and noticing that, 



grAU) = 9T,N-i{U) ■ (4) 
one obtains that the energy U fluctuates according to gamma distribution: 

9tAU) = Y^^if^Uf-' exp(-/3t/), (5) 

where (U) = TN and Var{U) = (U^) - {Uf = T'^N{N^1) - T'^N'^ = T'^N. 
• T = 1/ {3 and U are fixed and the multiphcity fluctuates. In this case 
one first writes a cumulative distribution function for the probability density 
function given by Eq. 



iV-l 1 

gtAu) = 1 - e j--r-xm~'^M-m- (7) 



* In the relativistic regime, where masses are neghgible in comparison with momenta, p « m, and 
for the cyhndrical phase space in which px^ >> px, we have E = p^. Therefore, when replacing discrete 
values Ei in the energy distribution function, Eq. ([2]), by a continuous variable E (and replacing sum- 
mations by integrals), we can put, for simplicity, the density of states p{E) = 1. In this case the energy 
distribution function ([2]) is normalized to unity. This simplification leads to simple interpretation of 
temperature parameter, T = (E), and specific heat, c = 1. 
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For energies {£'i=i,...,iv} such that 

N N+1 

Y.E,<U< E E,, (8) 
the corresponding multiphcity distribution (notice that U/T = {N)) is 

qtAN) = Gt,n+i{U) - GtAU) = ^ exp(-(iV)), (9) 



i.e., it has form of the Poisson distribution with 

Var{N) _ 1 



(10) 



• and total energy U are fixed and temperature T fiuctuates. Inverting dis 
tribution gT^N{U), Eq. (E]), one gets 



9uAT) = ^ j^dU'sTAU') = T^A^Uf-' exp(-/3[/), (11) 



013 Jo ^"'^^^ ^ r{N) 
Var{T) _ Var{(3) 1 



(12) 



(r)2 (/3)2 N' 

where = N/Ud.ndVar{P) = = N{N+l)/U'^-NyU^ = N/U'^. 

Notice that all hmiting distributions, Eqs. (Ej), (ED and (fTTI). have the formally 
identical form of a gamma distribution, 

UqtAU) = NgrAN) = pQuA^) = exp(-/3t/), (13) 



and also have identical respective relative fiuctuations, 

Var{U) Var{P) Var{N) 1 



{UY W {NY (N) 

For large Eq. ( fT3| ) becomes a Gaussian distribution. 



(14) 



Usually the temperature fluctuations in a system are related to its heat capacity under constant 
volume |8|9j . Cy, by = l/Cy and Cy = c{N) (in what follows we shall put specific heat c = 1). It 



means therefore that uij^ = uj^ = 1/{N) 
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which (with /i = N/9 and cr^ = N/9'^) is the distribution usually used to describe 
fluctuations in statistical physics [|8]. 



2.2 Generalized fluctuations 



Our considerations can be generalized by resorting to Tsallis statistics [5] with 
fluctuating U and T. In PJ7] it was shown that fluctuations of temperature in a 
heat bath in the form of a gamma distribution, result in a Tsallis distribution, 



hq{E) = exp 



q 



T 



E 
T 
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(16) 



with one new parameter, a nonextensivity parameter q {q < 2). The Boltzman 
distribution, Eq. (Ej), is recovered for g — > 1. It turns out [B] that the nonex- 
tensivity parameter q in Eq. ( |TB| ) is given by these fluctuations of temperature 
T: 



Var{T) 



= Uj^ = q — 1. 



(17) 



The further consequence of using Tsallis statistics is that now the joint N' 
particle Tsallis distribution with energies {£"4=1.. /^^}, 



})=c 



N 



Ejli E^ 
T 



1 
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(18) 



does not factorize into single particle distributions as g ({£'i=i,...,Ar}) in Eq, 
[TT]. As a result, the corresponding multiplicity distribution, which in the case 
of Boltzman-Gibbs statistics has a Poissonian form, cf. Eq. (Ej), now takes a 
Negative Binomial (NB) form [|TT], 



P{N) = 



T{N + k) 

r(iv + i)r(/c)" 



(19) 



where parameter k is given by the parameter q from Eq. (1161). 

1 



k = 



(20) 
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(in what follows we only consider the case of q > 1). On the other hand, from 
the definition of NB distribution 



l_VW{iV) 1 _ 1 

It means that fluctuations of and T are not independent, but related in the 
following way: 

^N-J^ = ^T- (22) 

However, NB multiplicity distribution can be obtained also as a result of fluc- 
tuations (in the form of a gamma function) of the mean multiplicity, N = {N), 
in the Poisson distribution. That is because in this case |12] (cf., also p!3lfT4] ): 

Pfm- rr^Af^^ 7^A^^-V^^ _ r{k + N) 7^ .... 
^^"io"^ iv! ■ r{k) "r(i + iv)r(/c)'(7 + i)^+^ ^ ^ 

that, for 7 = coincides with Eq. (HM. Therefore, in addition to Eq. ( 1211 ) 
one also has that 

1 _ Var{N) 

(iV)2 • ^^^^ 
For N = const we have a Poisson distribution [k = oo). Fluctuating 1/T accord- 
ing to gamma distribution and keeping U = const results in Var{N) / (N)'^ = uj^ 
and we recover Eq. ([22]) • Analogously, fluctuating U while keeping T = const 
gives us Var{N)/ (N)'^ = tufj. 

Fluctuating both U and T (and taking into account that N = U/T) one has 
that 



. Var n ■ 



2 



Var{U) Var{T) _ Cov{U, T) 



(25) 



{UY {TY {U){T) 
or, in terms of the scaled variances introduced before, 

9 9 9 / \ 

(jj^ = ujjj + LUj. — 2paj[/CJT, (26) 



where p = p{U, T) is the correlation coefficient {p G [—1, 1]) 



~k ~k •k 



***A similar relation connecting variables V , P and T, a;p = ojy +uj^, is known for almost a century 

m. 
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Comparing Eqs. ( flTh and ( HM and accounting for ([2SD one gets the following 
general relation between all fluctuating variables: 



(N) 



ufj — 2pujuujT- (27) 



This relation, which is our main result, generalizes Linhard's thermodynamic 
uncertainty relation given by Eq. ([T]). 

A word of explanation is in order. The use of | • • • | makes our formula ( 1271 ) 
general, i.e., valid for both cj^ > 1/ (N) and for cj^ = if = const. Actually, 
when all variables fluctuate one cannot have fluctuations of N smaller than the 
Poissonian. Observation of sub-Poissonian fluctuations, which would correspond 
to the case g < 1, always signal the presence of some additional constraints 
(like conservation of some quantum numbers, for example charges, cf., jl5]). 
We restrict ourselves to the case q > 1 and do not describe the region < 
uj% < 1/{N). We could, therefore, alternatively write Eq. ( I27| ) as 



- (1 - 2(5 (cj^j) =Uu + Ut- 2pujuUTj (28) 



where S{x) is Dirac delta. However, in what follows, we shall use Eq. (I27j) 



It is straightforward to see that when two of three variables are fixed results ob- 
tained using relation ( 1271 ) coincide with those obtained using Boltzman statistics 
(cf. Eq. (ED). When only one variable is kept constant we have: 



• for U = const one gets Eq. (I22D obtained using Tsalhs statistics; 

• for T = const one gets 

- ^/{N) = ujI, (29) 

i.e., fluctuations of energy result in fluctuations in multiplicity identical to 
those induced by fluctuations in temperature, cf. Eq. 



• for A = const one gets 1/A = + ojy — 2puJuUJT, where for = one 
obtains Eq. ([ID proposed in because — 1 < p < 1 the relative fluctuations 

2 2 

are constrained by the relation (cuu — lut ) < 1/A < {luu + ujt) . 
For the case when all variables are free to fluctuate we have Eq. ( HTh which can 
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be rewritten as: 



UJ 



N 



(N) 



= {uju — oJrf + 1uju^t{'^ — p)- 



(30) 



A Poissonian distribution of multiplicity (i.e., uj% = 1/{N)) is possible only for 



p = 1 when iou = lot- For p = one has 



N 



{NY 



UJ 



u 



+ Wy. Using Eqs. 



and (1211 ) one can express the nonextensivity parameter q by the respective 
fluctuations and correlations and write 

\q-l\ = {uju - urf + 2uuujt{1 - p) = co^ [(1 - + 2^(1 " P)] , (31) 

where ^ = lou/ut- Dependence of \q — l|/cjy on the correlation coefficient, 
/?, and relative fluctuations, ^, is shown in Fig. [H When U = const (i.e., ^ = 




1.0 1.5 



Fig. 1. (Color online) Dependence of the ratio \q — on the relative fluctuations, ^, for different 

values of correlation parameter p. 

0) we have the same situation we encountered for a fluctuating temperature, 
namely that \q — l\/ujj^ = 1. When we allow energy to fluctuate and add these 
fluctuations, \q — changes as shown in Fig. [H Notice that condit ion q < 2, 

implies that ojj/ + < 1, for 7^ 0, and that ojfj + oj^ < 1 for = 



~k ~k irk 



•k-kirk When energy U fluctuates, the pairs of variables, {U,N) and {U,T), cannot be indepen- 
dent simultaneously because Var{U) = {T)Cov{U, N) + {N)Cov{U,T) . This relation arises from a 
comparison of Eq. ()25p ) with the analogous formula evaluated for variable T = U/N. 
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3 An application 



As an example we compare fluctuations extracted from the distribution of dif- 
ferent observables in a high energy multiparticle production process. It should 
be remembered that Eq. (ETj) connects fluctuations of different observables, but 
defined in the same fragment of allowed phase space, whereas available data 
usually refer to different parts of this phase space. Therefore, corresponding q 
parameters are usually difficult to compare. For example, q = Ql obtained from 
rapidity {y) distributions, dN/dy, defined in so-called longitudinal phase space, 
are comparable with q evaluated from the multiplicity distributions, P{N), 
which are defined in the full phase space p!^6j . On the other hand, transverse 
momentum [px) distributions, dN/dpr, defined in the so-called transverse space, 
are described by much smaller values of q = qr- A first attempt to explore the 
relation ([T]) in high energy multiparticle production processes was presented in 

m- 

So far there are no data which would necessitate the use of nonzero correlations. 
In the case of uncorrelated fluctuations {p = ), one gets from Eq. ([23), using 
(Ell), that 

i = g - 1 = cj^ + . (32) 



In Fig. [2] we plot the energy dependencies of the nonextensivity parameter q ob- 
tained from different sources: from the multiplicity distributions, f{N) = P{N), 
i.e., from the full phase space |[IB] and from a different analysis of transverse 
momenta distributions, / (pr) (i.e., from the transverse phase space p!9ll20ll2T] ). 
The characteristic feature seen there is that, whereas the former show substan- 
tial energy dependence (and essentially follow results for q = qi obtained in 
[TBj from the analysis of dN/dy), the latter q = qr are only weakly dependent 



on the interaction energy (notice also that qris) from different sources plotted 
in Fig. [2] are roughly the same). To somehow compare q and qr, one needs 
some additional input. Assume that, being only weakly energy dependent, qr 
are also roughly independent of the energy fluctuations. It is now natural to 
expect that transverse characteristics are mainly governed by the fluctuations 
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1.40 n 

q from 
1.35- o f{N) 

1 30 • Pt) [Wibig] 
. ■ f(p,)[CMS] 
1.25- A f(^^)[NA49] 

1.20- 
1.15- 
1.10- 
1.05- 
1.00- 
0.95- 
0.90 -- 



or 




I — 

10 



-l 1 — I — I I I I 1 1 

100 



-I — I I I I I I 1 1 — I — I I I I I I 

1000 10000 



s'" [GeV] 



Fig. 2. (Color online) An example of energy dependencies of the nonextensivity parameters g obtained 
from different observables. Open symbols: q obtained from multiplicity distributions P{N) [181 (fitted 
by g = 1 + 1/k with 1/k = —0.04 + 0.0291ns, cf. [18j). Closed symbols: q = qr obtained from a 
different analysis of transverse momenta distributions, / (pt)- Data points are from, respectively, [19] 

for [Wibig], [20j for [CMS] and from [2T\ (data on fir = ^Jml + p\) for [NA49]). The dotted line 

represents a fit obtained in [T9| {qx = 1.25 — O.SSs"*^'*^^^) and the full line comes from our Eq. (I35p . 

of temperature, i.e., that we can write 

- 1 = ^ujI (33) 

(it is assumed here that fluctuations of temperature contribute equally to each 
of the components of momenta, hence the factor 2/3). Following this line of 
thought, i.e., assuming additionally that fluctuations of energy U are entirely 
given by its thermal part, one can write that in this case oj^ = and that 
both parameters, q and qx-, are connected by the following relation: 

g-l = 3(gT-l) or qT=^{2 + q). (34) 
A few explanatory remarks are in order. Namely, in p!9ll20] the pr distributions 
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were fitted using Tsallis distributions using the power 1/(1 — Qt)- However, one 
should keep in mind that in px distributions one has the pref actor px/T'^, not 
l/T, present in the usual exponential distributions. This fact results in a slightly 
different power, gr/(l — qr) [H] (this is the situation similar to the change from 



the so-called superstatistics A, in which only expression under the exponent is 
subjected to fluctuations, to superstatistics B, in which one fluctuates also the 
prefactor jTT], cf. also |22])- As a result the relation between qt and q is slightly 



modified and Eq. ( jMj ) now becomes 



1 + 2g 

Qt = ^7^- (35) 

2 + q 



Using values of q = 1 + 1/k = 0.896 + 0.0291ns obtained from P{N) 
the evaluated qr is shown in Fig. [2] (solid line) and compared to qr extracted 
from transverse momenta distributions /(pr) P^HIEUU^ . This, in turn, should 



be compared with the dotted line representing the fit in [19] resulting in qx = 
1.25 — 0.33s~^-^^^. Notice the good agreement of Eq. (ESj) with data which, in our 
opinion, justifies the statement that Eq. (ES) represents a kind of rule connecting 
fluctuations in different parts of phase space (modulo additional assumptions) . 



4 Summary 



Using nonextensive statistics applied to ensembles in which the energy (£"), tem- 
perature (T) and multiplicity (A^^) fluctuate, we have derived a specific relation 
connecting all fluctuating variables, Eq. ([23), whi ch genera 
modynamic uncertainty relation given by Eq. ([Tl) l^^^^^ 



izes Linhard's ther- 
This is illustrated 

using example taken from the multiparticle production processes. A possibility 
of connecting fluctuations appearing in different parts of phase space is indi- 
cated, cf. Eq. (ESj). 



•k-k-k-k-k Actually, the nonextensive approach used here is still subject to a debate about whether 
it is consistent with the equilibrium thermodynamics [23]. In this respect we would like to say that 
recently it was demonstrated on general grounds [U] that fluctuation phenomena can be incorporated 
into traditional presentation of thermodynamic and that the Tsallis distribution p6|) belongs to the 
class of general admissible distributions which satisfy thermodynamical consistency conditions and are 
a natural extension of the usual Boltzman-Gibbs canonical distribution Still other justification of 
a nonextensivity approach can be found in [25] . 
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